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Abstract 

y  In  some  spread  spectrum  environments,  the  low  energy  of  the  transmitted 
signals.  In  conjuctlon  with  the  existence  of  channel  noise,  do  not  allow 
the  distinction  between  collisions  and  lack  of  transmissions.  For  such 
environments,  and  for  the  Poisson  user  model,  we  propose  and  analyze  stable 
full  feedback  sensing  and  limited  feedback  sensing  synchronous  transmission 
algorithms.  -WCT  assume  binary  SNS  (success  versus  nonsuccess)  feedback  per 
slot,  and  the  possibility  of  transmission  of  phony  data  by  a  central  node. 
The  highest  throughput  attained  by  both  the  full  feedback  sensing  and  the 
limited  feedback  sensing  algorithms  is  0.322,  while  the  latter  induces  some¬ 
what  higher  delays.  This  is  compensated  by  the  robustness  of  the  limited 
feedback  sensing  algorithm  in  the  presence  of  feedback  errors  (in  contrast 
to  the  full  feedback  sensing  algorithm) ,  and  its  modest  requirements  on  the 
sensed  feedback  history.  j  .  7^* 
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1.  Introduction 

We  consider  multi  user  spread  spectrum  systems,  for  data  transmission.  In 
such  systems,  low  energy  data  packets  are  transmitted  to  a  central  node,  which 
has  broadcasting  capabilities.  Let  us  require  that  the  time  of  the  transmission 
channel  be  slotted,  where  each  slot  accommodates  a  single  packet.  The  users  of 
the  system  act  Independently;  thus,  simultaneous  transmissions  within  a  single 
slot  may  occur,  resulting  in  collisions.  In  addition,  due  to  the  low  energy 
of  the  data  packets,  in  conjunction  with  the  presence  of  noise  in  the  transmission 
channel,  make  collision  slots  indistinguishable  from  empty  slots  (slots  unoccupied 
by  transmissions).  That  is,  the  central  node  can  distinguish  only  between  success¬ 
ful  slots  (slots  occupied  by  a  single  packet)  and  nonsuccessful  slots.  The  result¬ 
ing  broadcast  feedback  per  slot  is  then  binary  SNS  (success  versus  nonsuccess) . 

We  consider  two  classes  of  multi  user  spread  spectrum  systems,  whose  common 
characteristics  are  as  described  in  the  above  paragraph: 

(i)  Systems  consisting  of  a  number  of  well  identified  static  users. 

(ii)  Systems  consisting  of  mobile  users.  In  the  first  class,  the  users  have  the  ca 
pability  to  maintain  in  memory  the  overall  broadcast  feedback  history.  In  the 
second  class,  and  due  to  their  mobility,  the  users  are  exposed  to  the  broadcast 
feedback  only  while  they  are  active;  that  is, from  the  time  they  generate  a  data 
packet,  to  the  time  that  this  packet  is  successfully  transmitted.  Therefore, 

full  feedback  sensing  transmission  algorithms  are  applicable  in  the  first  class, 
while  only  limited  sensing  such  algorithms  may  be  considered  for  the  second 
class.  For  both  classes,  we  will  initially  adopt  the  Poisson  user  model  (in¬ 
finitely  large  number  of  independent  Bernoulli  users),  since  this  model  represents 
a  limit,  and  since  algorithms  designed  for  it  are  robust  in  the  presence  of 
changing  user  structures.  For  this  user  model,  we  will  propose  and  analyze 
both  full  feedback  sensing  and  limited  feedback  sensing  transmission  algorithms. 
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in  Che  presence  of  binary  SNS  feedback. 

Given  some  transmission  algorithm,  its  throughput  is  defined  as  Che  maximum 

expected  traffic  Chat  the  algorithm  maintains,  with  finite  expected  delays.  For 

Che  Poisson  user  model,  the  throughput  Is  some  intensity,  X  packets  |sloC,  of  the 

Poisson  traffic  process.  The  algorithm  is  called  stable,  if  its  throughput,  6  ,  is 

* 

larger  than  zero,  and  the  interval  (0,6  )  is  then  Che  stability  region  of  Che 
algorithm.  The  delay  induced  by  the  algorithm  is  a  random  variable,  and  it  is 
defined  as  the  time  interval  between  the  instant  when  a  packet  is  generated,  to 
the  instant  when  it  is  successfully  transmitted.  The  expected  value  of  this 
variable  is  the  expected  per  packet  delay  induced  by  Che  algorithm.  In  packet 
networks,  the  algorithm  is  called  synchronous,  if  Che  channel  time  is  divided 
into  slots,  each  of  length  equal  Co  a  single  packet,  and  packet  transmissions  can 
only  start  at  the  beginning  of  some  slot. 

Let  us  consider  the  class  of  synchronous  stable  algorithms,  for  the  Poisson 
user  model.  Then,  among  the  existing  full  feedback  sensing  algorithms  within 
this  class,  Gallager's  algorithm  [Ij  provides  the  highest  throughput,  0.487,  and 
it  requires  ternary  (success  versus  collision  versus  empty)  feedback  per  slot. 

The  limited  feedback  sensing  algorithms  within  the  above  class  were  first  introduced 
by  Tsybakov  and  Vvedenskaya  [2].  In  [3],  higher  throughput  (0.42)  such  algorithms 
were  proposed  and  analyzed,  for  both  ternary  and  CNC  (collision  versus  noncollision) 
binary  feedbacks.  In  [5],  the  highest  throughput  known  limited  feedback  sensing 
algorithms  were  proposed  and  analyzed,  for  both  CNC  binary  and  ternary  feedbacks; 
the  throughput  for  CNC  binary  feedback  is  0.45,  and  the  throughput  for  ternary 
feedback  is  0.487.  In  contrast  to  the  algorithm  in  [1],  the  algorithms  in  [6]  are 
robust  in  the  presence  of  channel  errors,  and  they  induce  uniformly  good  delay 

characteristics.  In  [4],  a  unified  methodology  for  the  delay  analysis  of 
random  multiple  access  algorithms  is  presented.  This  methodology  has  been  used 
for  the  delay  analysis  of  the  algorithms  in  [3]  and  [5],  as  well  as  for  the  same 
analysis  of  Gallager's  algorithm  and  the  controlled  ALOHA  algorithm.  The  SNS 
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binary  feedback  in  the  full  feedback  sensing  environment  was  first  considered 
by  Berger  et  al  [6]. 

2.  The  Hodel-  A  Conmon  Algorithmic  Rule 

The  transmissions  from  the  users  are  assumed  synchronous.  That  Is,  each  user 
may  attempt  the  transmission  of  a  packet,  starting  only  at  the  beginning  of  some 
slot.  Time  Is  measured  la  slot  units,  and  slot  T  occupies  the  time  Interval  [T,T+1) . 
When  no  packet  Is  transmitted  In  slot  T,  the  slot  is  empty.  When  a  single  packet 
is  instead  transmitted  in  the  Interval  [T,T+1) ,  it  is  assumed  that  It  is  received 
correctly  by  the  central  node,  and  slot  T  is  then  a  success  slot.  Finally,  if 
at  least  two  packets  are  simultaneously  transmitted  in  [T,T+1),  a  collision  event 
occurs,  it  Is  then  assumed  that  the  Information  in  the  Involved  packets  is  lost, 
and  retransmission  is  then  necessary. 

It  is  assumed  that  the  central  node  can  distinguish  between  success  versus 
nonsuccess  (emptiness  or  collision)  events  per  slot,  and  that  it  broadcast  this 
SNS  binary  outcome  at  the  end  of  each  slot.  It  is  assumed  that  the  latter  broad¬ 
cast  Is  received  by  the  users  correctly,  and  without  delay.  Finally,  it  Is  assumed 
that  the  central  node  can  transmit  phony-noninformation  Including-packets,  whenever 
necessary.  Such  transmissions  are  necessary  to  help  users  distinguish  between 
collision  versus  emptiness  events,  with  some  delay.  Without  such  distinction, 
no  stable  transmission  algorithms  exist,  for  the  Poisson  user  model.  This  latter 
model  consists  of  an  overall  Poisson  traffic  process,  with  Intensity,  X  packets  [slot. 

In  the  sequel,  we  will  denote  by  x^,  the  broadcasted  by  the  central  node 
outcome  of  slot  T.  We  will  denote  by  x^=  S,  Che  success  event,  and  we  will 
denote  by  x^=NS,  Che  nonsuccess  event.  Also,  we  will  assume  that  the  central 
node  transmits  its  phony  packets  synchronously  (whenever  it  does) .  Whenever  such 
a  transmission  does  not  overlap  with  some  transmission  from  a  user,  a  success 
event  occurs,  and  the  corresponding  outcome,  S,  is  then  broadcasted.  Otherwise, 
a  collision  event  occurs,  and  the  outcome  NS  is  then  broadcasted.  Whenever  the 
transmission  of  a  phony  packet  by  the  central  node  is  required  by  the  operations 
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of  an  algorithm,  the  corresponding  algorithmic  rule  will  appear  as,  TRANSMIT  P. 

The  main  problem  In  the  presence  of  the  SNS  binary  feedback  Is  asynchronicity 
In  channel  history  Information,  among  the  active  users.  In  particular,  let 

■S;t£T-l,  and  "NS.  Then,  at  time  T+1,  the  active  users  are  divided  Into  two 
groups,  I  and  II.  Group  I,  If  nonempty,  contains  those  users  who  transmitted 
In  slot  T.  Group  II  contains  the  remaining  active  users.  At  T+l,  the  users  In 
group  I  know  that  T  Is  a  collision  slot;  the  users  In  group  II,  on  the  other  hand, 
can  only  deduce  that  T  Is  either  an  empty  or  a  collision  slot.  This  asynchronicity 
is  resolved  by  all  the  algorithms  in  this  paper,  via  the  followimg  common  actions 
in  slot  T+1. 

(1)  Each  user  who  transmitted  in  slot  T,  transmits  again. 

(11)  Each  user  who  did  ^ot  transmit  in  slot  T,  withholds  transmissions. 

(ill)  TRANSMIT  P. 

Due  to  the  above  actions,  it  is  clear  that  x^^j^“NS,  if  T  is  a  collision  slot,  and 
x^^j^=S,  if  T  Is  an  empty  slot.  Therefore,  at  time  T+2,  all  the  active  users  have 
the  same  information  about  the  outcome  of  slot  T,  at  the  expense  of  one  lost  slot. 

In  particular, 

1.  If  x^“NS  and  then  at  T+2  all  the  active  users  know  that  T  is 

a  collision  slot  (termed  x^=C) . 

2.  If  x^*NS  and  then  at  T+2  all  the  active  users  know  that  T  is 

an  empty  slot  (termed  x^=E) • 

If  T'^T+2  is  the  first  slot  after  slot  T+1,  such  that  x^,=NS,  then  the  actions 
(1),  (11), and  (ill)  are  repeated.  Thus,  the  following,  common  to  all  our  algorithms, 
algorithmic  rule  clearly  evolves: 


1  Rule  A 

i 

I 

i 

t 

i 


If  x^  j=NS,  for  all  j  such  that  05.j£2k,  where  k>0,  and  x^_2j^_j=S,then: 
(I)  In  slot  T+1: 

(a)  Each  user  who  transmitted  in  slot  T,  transmits  again. 

(b)  Each  active  user  who  did  not  transmit  in  slot  T,  withholds 
transmission. 

(c)  TRANSMIT  P. 


(a)  If  then  is  deduced. 

(b)  If  then  x^=E  is  deduced. 

(c)  Upon  deduction  of  the  outcome  x^,  the  outcome  is 

eliminated  from  the  memory  of  each  active  user.  The 
algorithm  operates  then  on  the  deduced  outcome  x^. 

Due  to  rule  A,  the  algorithms  in  this  paper  basically  operate  with  ternary 
feedback,  at  the  expense  of  some  slots  wasted  for  synchronization  among  all  the 
active  users.  Such  waste  clearly  results  in  throughput  reduction  and  delay  increase 
as  compared  to  the  corresponding  algorithms  which  operate  directly  with  ternary 
feedback. 

3.  Full  Feedback  Sensing  Algorithms-  The  Poisson  User  Model 

The  full  feedback  sensing  algorithms  require  that  each  user  know  the  channel 
feedback  history  at  every  point  in  time.  Due  to  that,  and  in  conjuctlon  with 
rule  A,  the  independent  users  can  then  act  synchronously,  to  resolve  collisions. 

In  this  section,  we  will  consider  two  algorithms,  for  the  model  in  section  2. 

The  first  algorithm,  called  FFSl,  is  Gallager's  algorithm,  as  explained  in  [7], 
in  conjuctlon  with  rule  A.  The  second  algorithm,  called  FFS2,  is  Capetanakis' 
dynamic  algorithm  [8],  in  conjuctlon  with  rule  A.  Both  algorithms  use  some  optimal 
arrival  interval  A,  and  in  their  analysis  we  will  denote  by  X  the  intensity  of  the 
Poisson  traffic  process. 


3.1  The  FFSl  Algorithm 

This  algorithm  (as  explained  in  [7]) utilizes  absorptions  on  the  arrival 

axis.  We  thus  define  the  parameters  L.,  and  W„  as  follows. 

N  N 

^N’  portion  of  the  initial  arrival  interval.  A,  that  is  resolved  by  the 

algorithm  before  absorption,  given  that  the  multiplicity  of  the  initial 
collision  in  A  is  N. 

^N"  expected  number  of  slots  needed  by  the  algorithm  to  resolve  the 

collision  in  the  portion  of  the  initial  arrival  interval  A,  before 
the  initial  collision  is  observed. 


From  Gallager's  algorithm  in  [71,  in  conjuction  with  rule  A,  in  section  2 


we  easily  conclude  that  the  following  recursions  hold. 


Lo  -  2  .  Lj  -  1 


4  ;  w.p.  0.5 


2  +  L2  ;  w.p.  0.5 


2  +  Ljj  ;  w.p.  2 


Ljj  -  <  3  +  ;  w.p.  N2 


2  +  L£  ;  w.p. 


W  ^  W  ^  X 

”0  1  ^ 


-1-1  /n\  -N 

2^*2  Vi  UP 


““  ■  -1 
;N>2  (  2  W. 


;  w.p.  (i)2"^  ,  if0,l 


From  (1)  and  (2),  we  respectively  obtain, 

^[1-2-N+1]-1(2+N2"\j^_^  +53  ;  N>2 


6  ;  N“2 


"n  ■  2'Ml-2*''|-ltN2-Vl  +  E  (l)^’\' 


;  N>2 


Let  us  define, 


X  =  XA 


L(x)  =  23  ^  e"’'  L 

k=0 


w(x)  =•  2Z  fr  e”"'  W 
k=0 
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Then,  the  stability  requirement  on  the  algorithm  clearly  Imposes  the  condition. 


L(x)  <  X  ^xW(x),  and  the  throughput,  X  ,  Is  then  computed  as. 


X  ■  sup[xW(x)L  ^(x)] 


If  X  Is  the  value  that  attains  the  supremum  In  (6),  then,  the  optimal  Initial 


*  *  -1 


*  ★ 


arrival  Interval,  A  ,  equals  x  (X  )  .  In  appendix  A,  we  compute  X  and  A  .  We  find. 


X  -  0.321946 


A  -  3.944 


We  will  perform  the  delay  analysis  of  the  algorithm  In  section  5,  together  with 


the  same  analysis  for  the  remaining  algorithms  In  this  paper. 


3.2  The  FFS2  Algorithm 


As  In  [8],  this  algorithm  also  utilizes  an  Initial  optimal  Interval,  A,  on  the 


arrival  axis.  Let  us  define. 


The  expected  number  of  slots  needed  by  the  algorithm  to  resolve  an  Initial  collision 
of  multiplicity  k,  before  this  collision  Is  observed. 


Then,  the  operations  of  the  algorithm  in  [8],  in  conjuction  with  rule  A,  In  section 


2,  clearly  induce  the  following  recursions. 


Lq  =  2  ,  =  1 


(2  +L 
U  L, 


;  w.p.  2 


From  (8),  we  obtain. 


2  ;  k-0 


1  ;  k-1 


[1_2~*^+^]'^  [2(1+2”’^)  +S(J)2"'^  L^I  ;  k  > 2 


Let  us  define. 


c)  »  e  ^ 


k!  H 


,  X  =  XA 


Then,  the  inequality,  L(x)  <  X  x,  represents  the  stability  condition  for  the 


algorithm,  and  induces  the  following  expressions  for  the  throughput,  X  ,  and  the  optimal 


initial  arrival  interval,  A  . 


8 


,*  *-.1*  * 

\  m  sup  Ix  L  (x)J  -  X  L  (x)  ,  &  -X  (A)  Cll) 

X 

*  * 

In  appendix  A,  we  exhibit  the  computation  of  the  parameters  A  and  A  .  We 
found. 


X*  -  0.30062 
A*  -  4.158 


(12) 


The  delay  analysis  of  the  algorithm  will  be  performed  in  section  5. 


4.  A  Limited  Feedback  Sensing  Algorithm  -  The  Poisson  User  Model 

The  limited  feedback  sensing  algorithms  require  that  each  user  follow  the 
channel  feedback  history  only  from  the  time  he  generates  a  packet,  to  the  time  that 
this  packet  is  successfully  transmitted.  That  is,  each  user  follows  the  feedback 
only  while  he  is  active.  In  this  section,  we  consider  an  algorithm,  named  LFS,  which 
is  basically  an  adaptation  of  the  LSTFA  algorithm  in  [5],  for  the  SNS  binary  feedback. 
This  adaptation  is  not  as  straightforward  as  with  the  algorithms  in  section  3.  The 
simple  superposition  of  rule  A,  in  section  2,  on  the  LSTFA,  applies  only  when  the 


latter  is  in  class  2  (see  [S]).  We  thus  present  here  the  full  description  of  the 
LFS  algorithm. 

4.1  The  Description  of  the  LFS  Algorithm 

The  algorithm  is  performed  by  each  user  independently,  and  at  each  point  in  time, 
it  distributes  the  newly  arrived  and  the  nontransmitted  packets  across  two  classes, 

A  and  5.  Transitions  in  time  within  and  across  the  classes  are  controlled  by  the  ope¬ 
rations  of  the  algorithm.  Class  A  contains  those  packets  which  can  not  yet  decide  if 
the  system  is  empty,  while  class  8  contains  those  packets  which  know  that  the  system  is 
empty.  All  packets  in  class  8  can  thus  act  synchronously,  while  the  packets  in  class 
A  can  not. 

Each  nontransmitted  packet  follows  the  rules  of  the  algorithm,  utilizing  a  set  of 
parameters,  R,A,L^,Tj^,T^,£,  and  G^.  Among  those,  parameters  R  and  A  are  system  para¬ 
meters,  and  they  are  subject  to  optimization,  for  the  satisfaction  of  the  desired  through¬ 
put  versus  expected  delay  tradeoff.  Parameters  L.,T, ,T  ,t,  and  G.  are  recursively 

A  1  ■  g  A 
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updated,  as  dictated  by  the  algorithmic  rules.  The  above  parameters  can  be  Inter¬ 


preted  as  follows. 

A  :  An  initial  arrival  Interval. 


R  :  An  upper  bound  to  the  number  of  consecutive  empty  slots  allowed  during  a 
collision  resolution  interval,  where  R  ^  1.  When  a  packet  observes  R  +  1 
consecutive  empty  slots,  it  knows  that  there  is  no  collision  resolution  in 
process.  The  pattern  corresponding  to  R  +  1  empty  slots  is  here  represented 
by  R  +  1  consecutive  slot  pairs,  each  with  outcomes  (NS,S). 

G^:  The  number  of  consecutive  empty  slots  observed  by  the  packet,  after  it  moves 
to  class  B.  The  consecutive  empty  slot  patterns  are  as  in  the  interpretation 
of  the  parameter  R. 


L^:  The  number  of  slots  containing  packets  from  class  A,  from  the  time  when  the 
packet  was  generated,  to  the  current  time.  If  the  slot  within  which  the 


packet  was  generated  contains  packets 

the  number  L. . 

A 

T  :  The  time  elapsed  from  the  time  instant 
^  current  time,  minus  the  examined  after 


from  class  A,  then  it  is  included  in 


when  the  packet  was  generated,  to  the 
the  above  time  instant  interval. 


t  :  The  total  length  of  the  arrival  interval,  which  is  transmitted  in  the 
current  slot. 


Tj^:  The  time  length  between  the  instant  when  the  packet  was  generated,  and  the 
ending  point  of  the  most  recent  arrival  interval  currently  chosen  for  trans¬ 
mission  (see  figure  1).  All  the  packets  in  the  interval  that  corresponds  to 
the  length  Tj^  belong  to  class  B. 

We  also  define  the  parameter  T,  as  follows. 


T  :  The  time  elapsed  from  the  time  instant  when  the  packet  was  generated,  to  the 
current  time. 


Upon  arrival,  each  packet  initiates  the  algorithm  independently,  following  the 
rules  below.  Simultaneously,  step  (I)  (c)  of  rule  A  is  followed  by  the  central  node, 
a.  Initialization 


Let  a  packet  arrive  at  the  time  instant  t  ,  t  e[T'-l,T').  The  packet  observes 

A  A 

then  the  feedback  x^,  and  continuously  observes  all  the  feedbacks  from  this  point  on 
until  it  is  successfully  transmitted.  At  T' ,  the  packet  moves  to  class  A  below,  with 
initial  values  T=T'-t  ,  L  =0,  and  G  =0. 


.0 


b.  Class  A; 

Each  user  has  a  register  with  2(R-fI)  spaces.  Each  user  with  a  packet  for  transmiss¬ 
ion,  who  belongs  to  class  A,  observes  the  feedback  time  sequence  and, 

1.  If  the  feedback  sequence  S,S,S  appears, then  he  sets  moves  to  class  8. 

2.  If  the  feedback  sequence  S,S,NS,S  appears,then  he  sets  moves  to  class  8. 

A 

3.  If  the  feedback  sequence  S,NS,NS,S,S  appears,  then  he  sets  moves  to  class  8. 

4.  If  R-fl  consecutive  feedback  pairs  (NS,S)  appear,  then  he  sets  L^b2R4-1,  and  moves  to  class  8. 

c.  Class  8; 

All  packets  in  class  8  act  as  follows: 

Start  with  T  “T,  G.“0,  and, 
g  ’  A  ’ 

2.  £  -  A,  Tj^  -  T-L^ 

Then, 

2.1  If  T  -(T-T,)<  ^ 

g  1  - 

2.1.  (a)  Set  T  •♦■T  +1,  set  T-»T+1,  set  G.  -0  and  TRANSMIT 
S  g  A 

2.1.1  If  X^**S,  the  packet  is  successfully  transmitted. 

2.1.2  If  X^»NS,  then  set  L^-*L^+2, 

2.1.2.  (a)  Set  ,  and, 

2. 1.2.1  If  Tg-(T-Tj^)^ -t,  move  to  step  2.1(a) 

2. 1.2. 2  If  T  -(T-TJ>£,  set  T -»-T  +1,  set  T-^T+1,  and, 

g  1  g  g 

2. 1.2. 2.1  If  X^-S,  set  Tg-*Tg-^,  set  L^-»L^+1,  and  move  to  step  2.1.  (a) 

2. 1.2. 2. 2  If  X,_-NS,  set  T  -►T  +1,  set  T-^T+1,  set  L^-»L^+1,  and, 
t  g  g 

2. 1.2. 2.2.1  If  X^“NS,  set  move  to  step  2. 2. 1.2. 

2. 1.2. 2. 2. 2  If  X„=S,  set  L.-^.+l,  set  T  -*-T  -  I, 

T  A  A  g  g 

set  and. 

If  G^<2R,  move  to  step  2.1.2. (a) 

If  G^-2R,  move  to  step  2.1. (a) 

2.2  If  T  -(T-T,)>£.,  set  T  -►T  +1,  T-*-T+l,  and, 

g  1-  g  g 

2.2.1  If  X  -NS,  set  T  -►T  +1,  T-^T+2,  L.-L.+l. 

T  g  g  A  A 


2. 2. 1.1  If  X_-S,  set  T  -*1  -I,  g  =G.+2,  and 

1  e  e  A  A 
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update  as  follows:  L^-*-L^+l.  Then, 

If  then  L^“3 

If  Chen  L^“3 

If  L^-2R+3,  then  L^-2R+1 
and  move  to  step  2. 

2. 2. 1.2  If  X.j,-NS,  set  +1,  set  ^^"0,  and, 

2.2. 1.2.  (a)  Set 

2. 2.1. 2.  (b)  Set  T  -►T  +1,  set  T>T+1,  and, 

S  S 

2. 2. 1.2.1  If  X_-S,  set  T  -*7  -I,  set  L.-»-L.+l 
I  R  R  A.  A 


set  T  +1,  set  T^T+1,  set  G.-O,  and, 
g  8  A  ’ 

2. 2. 1.2.1. 1  If  X_*»S,  set  T  ■♦T  ,  set  L,«4,  and  move  to  step  2. 

X  g  g  A 

2. 2. 1.2. 1.2  If  X_-NS,  set  T  -^T  +1,  T^T+l,  L,-^L,+1. 

X  g  g  A  A 

2. 2. 1.2. 1.3  Move  to  step  2. 2. 1.2 

2. 2- 1.2. 2  If  X_-NS,  set  T  -►T  +1.  T-^T+1,  L^-^^L.+l 
X  g  g  A  A 

2. 2. 1. 2. 2.1  If  X^»NS,  set  L^-^L^+1,  and  move  to  step  2. 2. 1.2. (a) 

2. 2. 1.2. 2. 2  If  X_*S,  set  T  -►T  — t,  set  L,-»-L,+l,  set  G  =»G,+2  and, 

T  ’  8  g  A  A  ’  A  A 

2. 2. 1.2. 2. 2.1  If  G,<2R,  move  to  2. 2-’. 2. (a) 

A 

2. 2. 1.2. 2. 2. 2  If  G  -2R,  move  to  2. 2.1. 2. (b) 

A 

2.2.2  If  X^“S,  set  *****  update  as  follows; 

If  L^-3,  then  L^-4. 

If  L^-2R+1,  then  L^-2R+2. 

Otherwise,  set 

Set  and, 

2. 2. 2.1  If  T  -(T-T, )<£,  move  to  2.1. (a) 

g  1  - 

2. 2. 2. 2  If  T  -(T-T,)>£,  set  T -.-T  +1,  T.4-T+1  and, 

g  1  g  g 

2. 2. 2. 2.1  If  X^=S,  set  and  move  to  step  2. 

2. 2. 2. 2. 2  If  X_-NS,  set  L,-*-L  +1,  T -♦T  +1,  T-«-T+l. 

X  A  A  g  g 

2. 2. 2. 2. 2.1  If  X^=NS,  set  h^L^+1  ,  and  move  to  step  2.2. 1.2. 

2. 2. 2. 2. 2. 2  If  X.j^“S,  set  s**'*  move  to  step  2. 


2 


In  figure  2,  we  include  the  flowchart  of  the  LFS  algorithm.  We  point  out  that  the 
algorithm  is  a  modification  of  the  FFSl  algorithm.  The  core  element  in  the  modification 
is  the  initial  observation  that  in  the  FFSl,  a  collision  resolution  interval  ends  with 
two  consecutive  successful  transmissions.  Based  on  this  observation,  a  packet  operating 
with  the  LFS,  recognizes  that  the  system  is  empty,  if  one  of  the  following  feedback 
sequences  occurs;  (i)  S,S,S  (ii)  S,S,NS,S  (ill)  S,NS,NS,S,S.  The  first  sequence  corresponds 
to  two  consecutive  successes.  The  second  sequence  corresponds  to  two  consecutive  successes, 
followed  by  an  empty  slot.  The  third  sequence  provides  the  information  that  two  consecutive 
successes  occured. 

4.2  System  Stability 

Let  us  consider  the  evolution  of  the  LFS  algorithm,  as  seen  by  an  outside  observer. 

Let  the  time  T  be  measured  in  slot  units,  and  let  the  operation  of  the  algorithm  start  at 
T=0.  At  T-2R+2,  there  will  be  then  2R+1  slots  containing  packets  from  class  A,  and  one 
slot  containing  packets  from  class  8.  Let  us  define  the  variables  T^,  D^,  L^(A) ,  and  I^, 
as  follows. 

Tg«2R+2 

T  :  The  first  time  after  T  , ,  such  that  the  distribution  of  the  packets  in  the 
n  n-l 

;n^l  unexamined  Interval  is  Poisson. 

D  :  The  total  length  of  arrival  intervals  containing  packets  from  class  B,  at  time 
n 

;n>0  T  ,  where  D=l. 

—  n  U 

L  (A) :  The  number  of  slots  containing  packets  from  class  A,  at  time  T  . 
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The  triple  S  "(D  ,L  (A),I  )  describes  the  state  of  the  system  at  time  T  ,  as  induced 
n  n  n  n  n 

by  the  operation  of  the  algorithm,  and  the  sequence  {S^}  is  a  Harkov  chain.  In  addition, 

2  <  L  (A)  <  2R  ;  Vn  >  0 
—  n  —  — 

D  >0  ;  Vn  >  0 

n  —  — 

and  the  values  of  D  are  demmerable.  Given  n,  given  a  state  value  s  *  (d  (A), I  ),  then, 

n  n  n  n  n 

^  A,  an  arrival  interval  of  length  A  is  chosen  by  the  algorithm  for  transmission. 

If,  on  the  other  hand,  d^  <  A,  then  no  packet  knows  the  value  d^,  and  each  packet  in  class  8 

selects  •  A.  Let  us  define, 

H  -  T  -  T  (13) 

n  n+1  n 

Then,  the  set  includes  1.1. d.  algorithmic  sessions,  and  as  in  section  5,  reference 

[5],  it  can  be  shown  that, 


E{H  |d  <  A)  <  E{h  Id  -  A}  -  e{H  |d  >  A}«»  ;V  A  <  “ 
n  n  “  n  n  n  n  ^ 

Consider  the  time  T  ,  let  the  value  of  the  system  state  be  then  s  •  (d  (A), I  ), 

n  n  n  n  n 

and  let  us  denote  A  ■  min(d  ,A).  At  time  T  . , ,  the  algorithm  examines  then  a  subset, 

n  n  u+i. 

^n’  ^n’  transmission.  For  as  in  (13),  and  following  the  same  procedure  as  in 

section  5  of  reference  (51,  we  conclude  that  for  stability  of  the  Markov  chain  it 

is  necessary  and  sufficient  that, 

EiagjsQ  -  (A,t(A),I)}  >  E{Hq1sq  -  (A,^(A),I)}  (14) 

Given  the  parameters  R  and  A,  the  throughput  X  (A,R)  of  the  LFS  algorithm  is  then 
the  maximum  intensity  of  the  Poisson  input  traffic,  that  maintains  the  condition  in 


4-^1 


(14).  Optimizing  then  with  respect  to  A,  we  conclude  that  given  R,  the  throughput, 

*  A 

X  (R),  of  the  algorithm  is  as  follows,  where  A  ^  denotes  the  optimal  initial  arrival 


interval. 


X  (R)  «  sup  X  (A,R)  ••  X  (A  _,R) 
A  ^ 


In  appendix  A,  we  exhibit  the  method  for  the  computation  of  the  throughput  X  (R) . 
In  table  1  below,  we  include  the  values  X  (R),  and  the  corresponding  values  A  ,  for 
various  values  of  the  parameter  R.  As  R  increases,  the  throughput  increases  as  well. 
The  increase  of  the  value  R  may  result  in  increase  of  the  expected  per  packet  delays. 


however.  This  last  issue  will  be  further  discussed  in  section  5,  where  the  delay 
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analysis  of  the  LFS  will  be  Included.  The  throughput  X  (6)  Is  the  same  with  the  througliput 


R 

X*(R) 

* 

^R 

1 

0.293376 

4.05623 

2 

0.315648 

3.960099 

3 

0.32043 

3.90099 

4 

0.32158 

3.9803028 

5 

0.321871 

3.976739 

6 

0.321943 

3.97585 

Table  1 

The  Throughput  of  the  LFS  Algorithm 
5.  Delay  Analysis  -  The  Poisson  User  Model 

In  this  section,  we  study  the  expected  per  packet  delays  Induced  by  all  the  three 
algorithms  in  sections  3  and  4.  The  delay  analyses  for  all  the  three  algorithms  are 
based  on  the  following  powerful  regeneration  theorem.  The  same  theorem  was  used  in  [4] , 
for  the  delay  analysis  of  the  Gallager  and  the  controlled  ALOHA  algorithms,  and  was  also 
used  in  [3]  and  [5],  for  the  delay  analysis  of  the  corresponding  algorithms. 

Theorem  A 

Let  the  discrete  time  process  {X  }  .  ,  be  regenerative  with  respect  to  the  renewal 

n  n>l 

process  ■  ^l+l~^i’  length  of  the  ith  regeneration  cycle, 

and  let  £(.)  be  some  nonnegative  real  valued  measurable  function.  Then, 

A  A 

If  C  -  E{Cj^}  <  <»,  and  S  •  <  **  ,  we  also  have, 

1  ^  —1  —1 

^lm  n  f(X^)  -  lim  n  E{^Ej^f(X^)}  -  SC  ,  w.p.l 

n->oo  n-»«> 

If  in  addition  to  the  finiteness  of  C  and  S,  the  distribution  of  Cj^  is  also  aperiodic, 
then  X^  converges  in  distribution  to  X^,  and, 

E{  f(Xj}  -  SC~^ 

We  point  out  that  the  process  called  regenerative  with  respect  to  the 

renewal  process  every  positive  integer  M  and  every  sequence  tj^,..,tj^,  such 

that  0  <  t,  <  ...  <  t„,  the  joint  distribution  of  . K  j-d  independent  of  i. 

-  1  -  -  .M’  V^i 


The  random  variables  and  respectively  called 

regeneration  points  and  regeneration  cycles.  Under  the  conditions  in  theorem  A,  the 
limiting  average  and  the  mean  of  the  limiting  distribution  of  exist,  coincide, 

and  are  finite.  In  addition,  their  conmon  value  is  then  given  in  terms  of  the  per  cycle 
quantities,  S  and  C. 

Referring  now  to  the  three  algorithms  in  this  paper,  we  will  define  as 

in  theorem  A,  the  delay  process  induced  by  the  corresponding  algorithm.  In  each  case, 

we  will  show  that  this  process  is  regenerative  with  respect  to  a  renewal  process 

with  easily  identified  regeneration  points  Then,  we  will  use  theorem  A,  to  establish 

the  existence  of  steady  state  delays,  and  to  compute  the  first  order  steady  state  moments 

(expected  delays),  by  appropriately  selecting  the  function  f(*).  For  each  one  of  the  three 

algorithms,  we  will  replace  the  process  theorem  A,  by  the  process 

where, 

V  :  The  delay  experienced  by  the  nth  arrived  packet.  That  is,  the  time  between 
'''  its  arrival  and  its  successful  transmission. 

D  -  E{P  }  ,  D  -  lim  D  (16) 

n  n  n-*o»  n 

5.1  The  Delays  of  the  FFSl  Algorithm 

Consider  the  operation  of  the  FFSl  algorithm.  Let  T  be  a  time  instant  that  corresponds 
to  the  beginning  of  slot  T,  and  it  is  such  that  all  packets  that  arrived  in  (0,  t^) ,  where 
t^  is  some  time  instant  such  that  t^<T,  have  been  successfully  transmitted,  and  there  is 
no  information  about  the  arrivals  in  tt^,T).  Then,  T  is  called  a  collision  resolution  in¬ 
stant  (CRI),  and  d^=T-t^  is  called  the  "lag  at  T."  At  every  CRI,  the  algorithm  restarts 
itself  by  selecting  an  initial  arrival  interval  =  min(A,d^).  Let  us  define  the  sequence 

{T  }  of  time  points, as  with  the  "0.487"  algorithm  in  [4].  Specifically,  let  ^^^=1,  dj^-2, 

^  ill 

and  let  T^^j,  be  the  first  instant  (corresponding  to  the  beginning  of  some  slot)  after  T^, 
such  that  d  -2.  Let  R,  ,  i>l,  denote  the  number  of  packets  successfully  transmitted 

4+1 

in  the  interval  Then,  t^l,  is  the  number  of  packets  successfully 

transmitted  in  the  interval  (T.,T  ],  and  the  sequence  is  ^  renewal  process. 
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since  ^  sequence  of  nonnegative  i.l.d.  random  variables.  In  addition,  the 

delay  process,  {V  }  ,  of  the  algorithm  is  regenerative  with  respect  to  the  process 

n  n>i 

^^i^i>l*  applying  theorem  A,  with  function  f(»)  the  identity  function,  we  have  that, 

if, 

C  -  E{Cj^}  <00  and  S  -  *  then, 

1-1 

D  -  S  C"^  .  w.p.l.  (17) 

;where  D  Is  as  in  (16).  It  thus  suffices  to  compute  the  expected  values  in  (17),  and 
prove  that  they  are  'finite,  where  we  also  have, 

C  -  XH  ,  H  -  (18) 

We  call  the  expected  value  H  in  (18),  mean  cycle  length,  and  we  call  the  expected 
value  S  in  (17),  mean  cumulative  delay.  In  appendix  B,  we  exhibit  a  method  for  the 
computation  of  upper  and  lower  bounds  on  Che  mean  cycle  length  and  the  mean  cumulative 
delay.  Those  bounds  are  then  used  for  the  derivation  of  upper  and  lower  bounds-respectively 
denoted  and  -  on  the  mean  packet  delay.  D,  in  (16) .  The  latter  bounds  are  exhibited 
in  table  2  below,  for  various  Poisson  intensities,  within  the  stability  region  of  the 
algorithm.  We  point  out  that  using  the  same  method,  tighter  bounds  can  be  devised,  with 
additional  computational  effort. 


X 

0.01 

2.06868 

2.082206 

0.05 

2.44738 

2.52379 

0.09 

2.815948 

3.120604 

0.13 

3.244687 

4.02865 

0.21 

5.28434 

7.34515 

0.25 

10.03156 

14.98144 

0.29 

22.57134 

48.518118 

0.32 

289.41968 

599.96109 

Table  2 


Mean  Packet  Delays  for  the  FFSl  Algorithm 


For  this  algorithm,  the  CRTs,  and  the  sequences  ^*'i^i>l’ 

are  defined  exactly  as  with  the  FFSl  algorithm  in  section  5.1,  and  so  are  then  the 

quantities,  C,  S,  D,  and  H,  in  (17)  and  (18),  and  the  function  f(*).  Again,  theorem 

A  applies, if  the  mean  cycle  length,  H,  and  the  mean  cumulative  delay,  S,  are  both  bounded. 
Then,  the  mean  packet  delay,  D,  equals  X  ^S,  where  X  is  the  Intensity  of  the  Poisson 
traffic  process.  In  appendix  B,  we  derive  upper  and  lower  bounds  on  the  quantities  H  and 
S.  Through  them,  we  derive  the  upper  and  lower  bounds,  u  and  D^,  on  the  mean  packet 

delay  D,  for  various  values  of  the  Poisson  Intensity  X,  within  the  stability  region  of 

the  algorithm.  We  Include  those  bounds  in  cable  3  below.  Again,  tighter  bounds  can  be 
found,  via  the  same  method,  and  with  additional  computational  effort. 


X 

D^* 

■■ 

2.0885448 

2.0885665 

2.5160811 

2.5241058 

0.09 

3.0976268 

3.1320109 

0.13 

3.935235 

4.0316111 

0.21 

7.7381153 

8.296022 

0.25 

14.049199 

15.641169 

0.29 

67.385942 

79.315479 

Table  3 

Mean  Packet  Delays  for  the  FFS2  Algorithm 


5.3  The  Delays  of  the  LFS  Algorithm 

The  sequences,  ^^i^i>l’  function  f(.),  and  the  quantities 

D,  H,  and  S,  are  again  defined  as  in  section  5.1.  The  delay  process  again 

regenerative  with  respect  to  the  renewal  process  theorem  A  applies.  The 

mean  cycle  length,  H,  and  the  mean  cumulative  delay,  S,  are  of  course  computed  differently 

here,  and  if  bounded,  the  mean  packet  delay,  D,  equals  X  ^S.  In  appendix  B,  we  exhibit 

u  Z. 

the  method  for  computing  upper  and  lower  bounds,  D  (R)  and  D  (R),  on  the  mean  packet 


delay  D,  for  every  value  of  the  algorithmic  parameter  R.  In  table  4  below,  we  include 


computed  values  of  the  bounds  d“(1)  and  D  (1)  -  that  correspond  to  algorithmic  parameter 
R  equal  to  one  -  for  various  values  of  the  Intensity,  X,  of  the  Poisson  Input  traffic. 

The  computation  of  the  bounds  D  (R)  and  0^(R),  for  R^  2  requires  Increasing  computational 
effort. 


X 

D^(l) 

d“(i) 

0.01 

3.05538 

3.2257 

0.05 

3.38136 

4.4241 

0.10 

3.9288 

6.6706 

0.15 

4.81452 

10.5137 

0.20 

6.6218 

18.4668 

0.25 

12.7931 

45.4634 

0.27 

23.8604 

94.0532 

0.28 

46.3069 

192.844 

Table  4 

Mean  Packet  Delays  for  the  LFS 
Algorithm  with  R*l. 

As  the  value  of  the  algorithmic  parameter  R  Increases,  the  mean  packet  delays  for 
small  Poisson  Intensities,  X,  increase  as  well.  Given  R^  and  R^  such  that  R^  >  R2, 
there  exists  some  Poisson  Intensity,  X°,  however,  such  that  for  Poisson  Intensities  above  X°, 
the  mean  packet  delays  Induced  when  the  logarithmic  parameter  Is  R^,  exceed  those  induced 
by  the  parameter  value  R^.  This  Is  shown  In  figure  3,  where  the  upper  bounds  d“(1)  and 
-  respectively  corresponding  to  R-l  and  R**2  -  are  plotted  against  X.  The  crossing 
point  is  Chen  approximately  equal  to  0.2.  That  is,  for  Poisson  rates  in  (0,0.2),  the  LFS 
with  R*1  is  superior,  while  for  such  rates  in  (0.2,  0.315),  the  LFS  with  R=2  is  recommended. 
5.4  Comparative  Discussion 

I 

In  figure  4,  we  plot  the  lower  bound,  D  ,  on  the  mean  packet  delay,  against  the 
Poisson  traffic  intensity  X,  for  the  FFSl,  FFS2,  and  LFS  with  R=*l  algorithms.  The 
upper  bounds  in  tables  3,4,  and  5  exhibit  similar  behavior.  The  FFSl  algorithm  induces 
the  lowest  mean  packet  delays,  for  every  X  value  within  its  stability  region.  It  is 


sensitive  to  channel  errors  (Leading  to  deadlocks),  however,  and  it  requires  full 
feedback  sensing.  The  LFS  algorithm,  on  the  other  hand,  is  robust  in  the  presence  of 
channel  errors  (as  the  algorithm  in  [5]),  and  only  requires  limited  feedback  sensing, 
at  the  expense  of  somewhat  Increased  mean  packet  delays.  Since  full  feedback  sensing 
is  impossible  in  mobile  spread  spectrum  environments,  the  latter  algorithm  is  then 
Indlspenslble. 


a.l 


A.l  The  IhrouKhput  of  the  FFSl  Aleorlthm 


We  use  the  following  lemma,  for  the  set  iL^jln  (1)  and  (3). 

Lemma  A1 

The  inequality  <  6k  holds,  for  every  k  ^  1. 

Proof 

We  prove  the  lemma  by  induction.  We  clearly  have  <  6k,  for  k«l,  2.  Accepting 
Chen  Chat  ^  6t  ;  and  using  (3),  we  easily  find  that  £  6(k+l) .  The 

proof  is  now  complete. 


Let  now  k  be  some  positive  integer,  and  let  us  then  adopt  Che  following  simple 
o 


bounds . 


1  <  L.  <  6k  ;  Vk  >  k 

—  V  —  —  r 


0<W,  <1  ;  Vk>k 

—  k  —  —  o 


(A.l) 


Using  (A.l)  and  (5),  we  then  obtain. 


e*B^(x)  -  X  52  W.  <  X  5!^  W.  -  e’^  x  W(x)  < 
k»0  ^  k«0 


w  It 


/‘^o  k  «  k  \ 

Vk-o*^'  k-k +i'"7 
k  ^ 

(«*  +  E  %  <Wk-l))  ^eV(,) 


(A.  2) 


e'‘L^(x)  -  i 


ok  *  k 

^  ^k  -  ^  TT  ^k 


ok  “»  k 

L(x)  1 E  V  ^k  ^  V 

k=0  k=k  +1 

o 


o  k 


6X6*  +  L  +  51  (L. -61^)  =  l'*(x) 

O  •  «  Ik  •  K 

k*l 


(A. 3) 


Then, 


B^(x)  ^  xW(x)  ^  B^(x) 
L“(x)  L(x)  L  (x) 


(A.4) 


'  .  ..n  O 

ScSR^^Sgk 

ri"! 

m 

1 

•  w  »»•  •• 

selected  k  -20,  we  computed  precisely  from  (3)  the  values  L,  ;  k  <  20,  and  we 

O  IC 


^  U  11  ^ 

then  found  numerically  the  suprema  of  the  ratios,  B^(x)/L  (x)  and  B  (x)/L  (x) .  Those 


suprema  were  both  attained  at  x  -1.266,  and  were  Identical  to  each  other  to  the  fifth 


decimal  point,  and  equal  to  X  In  (7). 


A. 2  The  Throughput  of  the  FFS2  Algorithm 


It  Is  easily  proven  that  lemma  A1  holds.  Then,  for  some  positive  Integer  k^,  and 


using  (10),  we  conclude  that  expression  (A. 3)  holds  here  as  well,  for  the  set  (L.) 


In  (9).  We  selected  k  -20,  and  we  computed  the  values  L,  ;  k  <  20  precisely  from  (9) 

o  K  — 


Then,  we  found  numerically  the  suprema  of  the  expressions,  x/V^ix)  and  x/L  (x) ,  where 


x/l'^Cx)  ^  x/L(x)  x/L  (x).  The  latter  suprema  are  both  attained  at  x  -1.25,  and  they 


are  Identical  to  the  fifth  decimal  point,  and  equal  to  X  In  (12) 


A. 3  The  Throughput  of  the  LFS  Algorithm 


Let  us  define, 


Q  -  A'^6, 


L.  -  E{H-  I  D  -  A,  k  packets  in  A) 


(A.5) 


-  E{n  I  packets  in  A} 


pk  A  (kj  2-k 


Then,  from  the  operation  of  the  algorithm,  we  deduce  the  following  recursions. 


I-O  ■  2.  Lj  .  1 


2+Zt  +  1  +  ;  w.p.  (Pq)^Pj  »  0  1  ^  1  R-1 


2+21  +  L. 


;  w-P-  °  1  ^  1  1  1  2 


(A. 6) 


2+2R  +  L. 


;  w-P-  (Pq) 


Wo  -  Wf  -  1 


23  2"J  +  a2"^"^  +  2"^"V  ,  ;  w.p.  (pj)  V.  0  <  f  <  R-1 

j=l  0  1  -  - 


^  2"J  + 

j-1  ^ 

23  2'J  +  2"*^  W 

J-1 


;  w.p.  (Po)^Pj  .  0  <  I  <  R-l.k  >  i  >  2 


,  k.R 

;  w.p,  (p^) 


(A. 7) 


From  (A. 6)  and  (A. 7),  we  respectively  find,  where  x  *0,  if  k<£, 

i-£  ^ 


;  k-0 
;  k-1 


.  k.R  k  I  f  .  ,  r  ^  k  k  I  T  1  ^  rii  n^ 

{l-(P(j)  “Pq  ~  T"  ■  )  \  2  S  Pi  ^i]  J  (A.8) 


1  ;  k-0,1 

|l-2“^pJ)’^-[2-pJ]'^l-2“^(pJ)®lp^|'^C2-pJr^[l-2’^(PQ)’^l. 

•  lP0^N^ 


Let  us  define. 


(A.9) 


L(x)  ^  X)  Jr  L 

k^O 

W(x)  »  X)  ®”’'  fr  L. 
k«0 


^  ;  where  0  ^  Wj^  j<  1;  V  k,  and  »diere  lemma  A.l  is  found  to  hold  on  the  sequence 

^  Then,  the  inequalities  in  (A. 2)  and  (A. 3)  hold  here  as  well,  and  so  does  (A.A) , 

■  where  again, 

S?  ,*  ,x  W(x)x 

I  *■'">“  ‘"lKT' 

^  For  all  the  R  values  in  table  1 ,  and  for  k^  -20  in  (A. 2)  and  (A. 3) ,  the  supreme 

9  of  the  bounds  in  (A.4)  where  indentical  Co  each  ocher  to  the  fifth  decimal,  and  equal 

A  .  A 

Co  Che  corresponding  values  A  (R)  in  cable  1.  The  values  were  found  as, 

KX  A  A  A 
r!  \-xrCA(r))/ 


b.l 


Appendix  B 


B.l  Delays  for  the  FFSl  Aleorithm 


Consider  the  CRIs,  T,  the  instants  t^»  and  the  Intervals  U^,  in  section  5.1. 

Then,  in  slot  T,  some  arrival  interval  [t^, »  is  transmitted,  where  after  some 

random  number,  of  slots,  another  CRI,  T'  is  reached,  with  a  corresponding. 

Let  us  define  the  following  quantities. 

:  The  number  of  slots  needed  to  resolve  an  arrival  interval  between  two 
subsequent  CRIs,  given  that  the  initial  length  of  this  interval  is  d. 

w^  :  The  length  of  the  actually  "examined*'  interval  before  absorption,  given  that 
the  length  of  the  InltlAl  Interval  is  d. 

P(‘£,w|d)  :  The  probability  that,  given  an  initial  length  d  Interval  between  two 

subsequent  CRTs,  the  actually  ''examined*'  before  absorption  interval  has 
length  w,  and  t  slots  are  needed  for  its  resolution. 


:  The  expected  number  of  slots  needed  to  reach  lag  equal  to  2,  starting  from 
lag  equal  to  d. 

The  operations  of  the  algorithm  induce  then  the  following  recursions. 


2  ;  w.p.  P(2,d|d) 

E{^d>  +  ’  ''•P* 


,  for  d  £  A 


(B.l) 


^d  "  **d-w*Z  •  ''****  P(^.w|A),  for  d  >  a 

Taking  expectations  in  (B.l),  we  obtain, 

{eUj}  +  Z  P(^.wld)  ;  d  <  A 

'“i’  *  ?  j  >  A 

4,w 


(B.2) 


The  expected  value,  H,  in  (18)  .  qual  to  H^.  where  is  eiven  bv  (B.2).  We 
thus  need  to  find  bounds  on  H2.  Following  the  same  procedure  as  with  the  ''0.487" 
algorithm  in  [4],  on  the  linear  system  in  (B.2),  we  find  the  following  bounds. 


^  E{^2)  +  A[2-2P(2.2|2)1  -P(l,2|2)  -  E{w2-42>1  c^(l-P(2.2|2)  -P(l,2|2)l 

+  P(1,2|2)(b  +c^{1-P(2,1|1)}1  £  H2  -  H  <  Ed^)  + 


+  A(2-2P(2,2|2)  -P(1,2|2)  -E{w2-£2>1  +c“(1-P(2, 2| 2)  -P'l,2|2)]  + 
+  P(l,2|2)rB  +c“{l-P(2,l|l)))S  H^* 


(B.3) 


b.2 


where, 

A  -  E{£^}  (E{w^}  -  eU^}]"^ 

p(d)  -  P"^(2,d|d)  {eUj}  +  ACeU^}  -  E{wj}l  -2AP(2,d|d)  | 

=  inf  p(d)  (B.4) 

l<d£A 

“  max  (-A,  sup  p(d)) 

l<d<A 

B  -  EUj^}  +  A[E{£j^}  -  E{wj^}  +  l-2P(2,lll)] 

From  (B.3)  and  (B.4),  we  computed  Che  upper  and  lower  bounds,  and  on  H, 
shown  in  Cable  B.l  below. 


X 

h“- 

0.01 

2.01435 

2.01731 

0.05 

2.12193 

2.14923 

0.09 

2.30954 

2.364348 

0.13 

2.514186 

2.671167 

0.21 

3.109886 

3.98608 

0.25 

4.852168 

6.60396 

0.29 

8.09345 

10.923344 

0.32 

111.98070 

161.02000 

Table  B.l 

Bounds  on  H  for  Che  FFSl  Algorithm 


Let  us  define  Che  following  quantities,  referring  to  the  beginning  of  this  appendix. 

N  :  The  number  of  packets  in  (t^,  t^,) 

ill  :  The  sum  of  the  delays  of  the  above  N  packets,  after  the  CRI,T. 

z  :  The  sum  of  the  delays  of  the  above  N  packets,  until  the  instant 

-  E{2|N-k,U^-l}  -  d"^  E{z|N-k,U^  -d}  ;  1  <  d  <  A 

-  E{i|»|N-k,U.j.-l}  -  E{(l;|N-k,U.j.-d)  ;  1  1  d  <  A 

Sj  :  The  expected  sum  of  delays  of  all  packets  transmitted,  when  starting  from 
lag  d,  the  algorithm  reaches  lag  2. 


b.3 

The  algorithm  Induces  then  the  following  recursions. 
0 


"d  -  ! 


•  t  *  1 


,for  d  £  A 

'  E{i|»|d}  +  E{z|d}  +  ;  w.p.  P(£,w|A)  ,  2 

-  E{i((|d}  +  E{z|d}  +  (d-A)N  +  ;  w.p.  P(£,w|A)  ,  for  d  >  A 

From  (B.S),  we  obtain. 


j  E{i(»ld}  +  E{zld}  +  £  P(-t.w|d)  ;  d  <  A 

)  ^.w 

m 

E{i|»|A}+  E{z|A}  +  (d-A)E{NjA}  +  Z  ;  d  >  A 

£,w 


1 


(B.6) 


;  where  for  S  as  in  (17),  we  have,  S  ■  82*  Following  the  same  procedure  as  with 


the 


u 

bounds  H  and  H  on  H,  we  find, 
<  $2  -  S  <  S“ 


;  tdiere. 


(B.7) 

.1. 


S“  -  G  +  c“  ri-P(2,2|2)]  +  P(1.2|2)  (s“-Mj^-M2-c^] 

S"^  -  G  +  c'^Il-P(2,2|2)l  +  P(1,2|2)Is^-M2-M2-c‘‘] 

G  -  E{<|)|2}  +  E{z|2}  +  2(Mj^+2M2)tl-P(2,2|2)l  +  M2E{(w2-^2)^>  "  (AM2-mj^)E{w2-^2> 
s®-  E{)|/|1}  +  E{z|l}  +  M2(i  +  ECwJ)  +  eUJ}  -  2EU2W^}  -  2[E{w2}  -  eUj^}] 

-  4P(2,l|l)J  +  M^(l  -  E{w^}  +  eU^}  -  2P(2,l|l)  |  +  c®[l-P(2,l|  1)  ]  ;  a-u,^ 

M2“  e{n|A}  2"^  tE{w^}  -  eU^}]"^ 

”1“  +  E{z|A}  +  M2lE{w^}  +  E{4^}  -  2E{w^4^} 


-  ae{n|a}| 

p(d)  -  P"^(2,d(d)  |E{4»|d}  +  M2E{(w^-4^)^}  -  (2M2d-Hl^)  E{w^-^^} 


(B.8) 


-  (4M2+2M^)  P(2,d|d)| 


o'*  «  sup  p(d) 
l<d<A 

I 

c  ■  inf  p(d) 
l<d<A 


^  u 

From  (B.S),  we  computed  the  bounds  S  and  S  ,  that  we  include,  for  various  values 


of  X,  in  table  B.2  below. 


X 

s'" 

0.01 

0.04173 

0.04194 

0.05 

0.26299 

0.26776 

0.09 

0.59920 

0.64864 

0.13 

1.12672 

1.31674 

0.21 

4.42339 

4.79694 

0.25 

16.5620 

18.17311 

0.29 

71.5080 

113.8769 

0.32 

14912.7540 

21498.9000 

Table  B.2 

Bounds  on  S  for  the  FFSl  Algorithm 

From  (17)  and  (18) ,  we  conclude  that  the  bounds  D**  and  D  ,  on  the  mean  packet  delay 
D,  are  computed  from  tables  B.l  and  B.2,  as  follows: 

d“  -  S“  (H^)"^  ,  -  X~^  (H**)"^  (B.9) 

B.2  Delays  for  the  FFS2  Algorithm 

Let  us  define  the  following  quantities. 

:  The  number  of  slots  needed  to  resolve  an  initial  arrival  interval  of 
length  d,  where  1  ^  d  A. 

P(£(d)  :  The  probability  that  ■  t. 

:  The  expected  number  of  slots  needed  to  reach  lag  2,  starting  from  lag  d, 
subject  to  successful  transmission  of  every  packet  in  the  arrival  interval 
of  length  min(d,A). 

S.  :  The  expected  sum  of  delays  of  all  packets  transmitted,  when  starting  from 
lag  d,  the  algorithm  reaches  lag  2,  and  all  the  packets  in  the  arrival 
interval  of  length  fflin(d,A)  are  successfully  transmitted. 

Let  the  quantities  N,  ij;,  and  z  be  as  in  section  B.l.  Then,  the  algorithm  induces 
the  following  recursions. 


“d*  I 


2  ;  w.p.  P(2|d) 


EU^  }  +  ;  w.p.  P(£|d)  ,1^2 

^d-A+£ 


,  for  d  <  A 


(B.IO) 


E{£^}  +  H,  ;  w.p.  P(£lA)  ,  for  d  >  A 


"4  «'>  W>  Ji"^ 


E{i(»|n-0}  +  E{z|N-0}  ;  w.p.  P(2|d) 


,  for  d  <  A 


^  *E{i|»ld}  +  E{zld}  +  ;  w.p.  P(^|d)  (B.ll) 

Sj-  E{il»jA}  +  E{z1a>  +  (d-A)  E{NjA}  +  ;  w-P-  P(^|A)  ,  for  d  >  A 

From  (B.IO)  and  (B.ll),  we  respectively  obtain. 


**d  I 


*d"l 


eU.)  +  Z  H,  P(-t|d)  ;  d  <  A 

EU^}  +  Z  PU\A)  :  d  >  A 

E{i|;|d}  +  E{z|d}  +  E  S.  P(4|d)  ;  d  <  A 

m  ^ 

E{ij>|A}  +  E{z1a}  +  (d-A)  e{n|A}  +  E  P(^1A)  ;  d  >  A 


(B.12) 


(B.13) 


Considering  the  quantities  H  and  S,  in  section  5.2,  we  have,  H  ■  H2  end  S  ~  82* 

Vie  are  thus  seeking  upper  and  lower  bounds  on  the  quantities  H2  and  82*  given  respectively 
by  (B.12)  and  (B.13).  Vie  note  that  here. 


;  where. 


P(4ld)  -  E  P(tlk)  P(kld) 
k 


P(k|d)  -  e 


P(Hk)  -  2'*"  lP(4-2lk)  +P(t-4|k)l  +  E  (h  2"'^  E  P(tJi)P(£_|k-i) 

1-1  ^  ^,+e,«^-2  ^  ^ 


k-1 


'1  2 


(B.14) 


To  compute  bounds  on  H.  and  8.,  we  use  iterations.  Let  and  8j^^  denote  the 

A  z  a  Cl 

corresponding  1th  iterations.  Then,  we  easily  conclude, 

.(i)  .  r  r„(l)  o(0) 


(B.15) 


»d  -  «d-A^£J  ^ 

t-1 


Md  +  <  Md  +  c“ 

—  a  — 


(B.16) 


b.6 


;  where, 

M  -  [A  -  E{£^} 

■  sup  f(d)  ,  «  inf  f(d)  (B.17) 

l£d<A  11A<A 

f(d)  -  P”^(2|d)  +  M[E{£^}-A]  -  2MP(2|d)j 

Using  K-40,  100  iterations,  and  via  (6.15),  (B.16),  and  (B.17),  we  found  the  lower 


and  upper  bounds,  H  and  h“,  on  H2  (and  thus  on  H) ,  shown  in  table  B.3  below. 


\ 

h“ 

0.01 

2.02178 

2.02178 

0.05 

2.14534 

2.14534 

0.09 

2.33792 

2.33793 

0.13 

2.63508 

2.63520 

0.21 

4.04580 

4.05625 

0.25 

6.40808 

6.47551 

0.29 

26.37999 

27.20669 

Table  B.3 


Bounds  on  H  for  the  FFS2  Algorithm 

For  the  quantities  S.,  we  started  with, 

a 

+  M,d  +  M„d^  <  <  M,d^  +  M,d  +  e“  (B.18) 

;  where, 

“  2'^  (A  -  eU^}]~^  E{n|A} 

-  [A-  Ea^}]"^  ^E{i|;(A}  +  E{z|A>  +  M2(A^  +  eU^}  -  ZAEU^}]-  Ae{n|A} 

e^  -  inf  f  (d)  ,  e“  -  sup  f  (d) 

l<d<A  °  l<d<A 

f^(d)  =  P~^(2ld)  ^E{t{)ld}  +  E{zld}  +  M2[eU^}  -  4P(2ld)]  + 

+  Mj^IeR^}  -  d  -  2P(2|d)]j 


(B.19) 


;  and  where , 

E{N|d}  -  Id  ;  lld<A 

E{l|»ld}  -  E  E{t;»lk}P(k|d) 
k 


E{z|d}  =  2~^ld^  ;  l<d<A 


E{i^lk}  -  - ^Zk-t-  Z  (j)2  ’‘(2E{i|;(i}+  (k-i)  E{f^  |  i} 


I.’V 


.-i.  j'-.' 


b./ 


E{l|»|0}  -  0  ,  -  1 


eU^}  -  2:  EU^Ik}  P(k|d) 
k 


E{£^|k)  -  [a  +  12. 2"*^  +  12.2"'^  EU^jk}  +  2  Z  (^)2'*‘[E^{4^|ra}  + 

V  m*l 

+  AEU^jm}  +  E{-d^|m}  EU^|k-m}]j  (B.20) 

The  iterations  for  S.  evolve  as  those  in  (B.15).  After  100  iterations,  we  found 

d 

Che  upper  and  lower  bounds,  S  and  S  ,  on  (end  thus  on  S) ,  Chat  are  Included  in 
Cable  B.A  below. 


X 

s“ 

0.01 

0.0A222 

O.OA222 

0.05 

0.26989 

0.27075 

0.09 

0.65178 

0.65901 

0.13 

1.3A811 

1.38107 

0.21 

6.591A3 

7.0A8A5 

0.25 

22.7A395 

25.057A9 

0.29 

531.6711A 

606.77905 

Table  B.A 

Bounds  on  S  for  Che  FFS2  AlgorlChm 

The  bounds  d“  and  are  compuCed  from  Cables  B.3  and  B.A,  and  from  Che  expressions 
in  (B.9). 

B.3  Delays  for  Che  LFS  AlgoriChm 

Here,  we  will  derive  bounds  on  Che  mean  packec  delay,  D,  when  Che  algoriChmic 
parameCer  R  equals  one.  The  meChodology  for  R^2  is  similar.  LeC  Che  quanCicies  and 

Sj  be  defined  as  in  seccion  B.l.  LeC  6q  and  d^  be  as  in  seccion  A. 2,  and  lec  be 
defined  as  in  secCion  B.l.  LeC  us  also  define,  where  T  is  some  CRI, 

N  ;  The  number  of  packers  in  6^. 

Nq:  The  number  of  packers  in  A  -  6^. 

P(£,wld)  :  The  probabilicy  chac  I,  and  chac  given  d=min(dQ,A),  Chen  6q»w. 
ijj  :  The  sum  of  Che  delays  of  Che  N  packeCs,  afCer  T. 

z  :  The  sum  of  Che  delays  of  Che  N  packets,  before  T-1. 


b.8 


Then,  Che  algorithm  Induces  Che  following  recursions. 


"d-vH-  I 


;  1  <  d  <  A 


(B.21) 


+  E  ;  d  >  A 

l,w 


E{z|d}  +  E{\|)ld}  +  E{N|d}  +  ECNpAQld)  +  E  P(^.w|d)S^_^^  ;  1  £  d  <  A 

'd  ^  "  (B.22) 

E{z|A}  +  e{(|>1a}  +  e{n1A}  +  X(d-A)  E{6q1A}  +  ECN^d^jA}  + 

+  E  ;  d  >  4 

t,v 

The  expressions  (B.21)  and  (B.22)  both  determine  infinite  dimensionality  linear 
systems.  Using  Che  methodology  for  the  delay  analysis  of  Che  algorithm  in  [5],  we 
conclude  chat  Che  systems  in  (B.21)  and  (B.22)  have  both  unique  solutions,  within  the 
class  of  quadratically  bounded  sequences,  if  E{6Q|dQ»*A}  >  E{.t^K  The  latter  inequality 
determines  the  stability  region  of  the  algorithm.  In  addition,  we  have  again  that  H=H2 
and  S-Sj,  where  H  and  S  are  the  quantities  in  section  5.3.  We  derive  upper  and  lower  bounds 
on  $2  and  (and  thus  on  S  and  H) ,  as  follows. 


<  H2  “  H  <  h“ 

<  $2  =  H  < 


(B.23) 


;  where. 


h"  -  E{£2>  +  2[1-P(2,2|2)]  f  +  c“[l-P(2,2(2)] 


/  »  h“  -  [c"*  -c^](l-P(2,2|2)] 


Mq  -  (E{5y|A}  -  E{£^}1~^  eU^} 


c^  =  max  (-M-  ,  sup  h(d))  ,  c^  =  inf  h(d) 

l<d<A  l<d<A 


h(d)  -  p"^(2,dld)|E{^^}  +  '  EiS^ld}]  -  2MpP(2,dld)} 

S‘'  -  E{z|2}  +  E{il<|2}  +  E{n|2}  +  E{Nq6q|2}  +  M2[A-4e~^^+  E{ (5q-£2) ^  1 2)  - 


-  4  E{6q-£2(2}]  +  M^[2-2e‘^^-  E{6q-£212}1  + 


u,,  -2X, 


% 


ct  -^1  /I  ~2Xs 

S  “  S  -  (e  -e  J  (1-e  ) 

M2  -  (E{6q|A}  -eU^}]’^X  E{6q|A} 

-  [E{6q|A}  -  EU^}1“^|e{z|A}  +  E{(j/|A}  +  E{nIA}  +  e{Nq6q|A}  - 

-  X  A  E{6q(A}  +  M2  E{(6q-£^)^1a}} 

=  sup  f(d)  ,  =  inf  f(d) 

l<d<A 


f(d)  =  p"^(2,dld){E{zId}  +  E{.()|d}  +  E{N|d}  +  E{Np6Q|d}  + 
+  M2[E{(5Q-£^)^(d}  -  2d  E{6g-£^|d}  -  4P(2,d|d)] 

-  M^[E{6Q-£^ld}  +  2  P(2,dld)][ 


(B.24) 


As  compared  to  the  computations  In  sections  B.l  and  B.2«  the  difficulty  here  is 

the  computation  on  the  expectations  included  in  the  expressions  in  (B.24).  Let  us  define, 

E^  {Xjd.k}  ;  The  conditional  expectation  of  the  random  variable  X,  given  that 

A  equals  A,  given  d,  where  d  £  A,  and  given  that  there  are  k  packets 
in  d. 

Then, 

E^  {X|d}  =  1]  E^{x|d,k}  e"^'^  (B.25) 

k=0 

Let  us  define  the  subset  V  of  (0,1],  as  follows, 

V  =  tv  :  v=l,  orv  =  2  ^  +  2  ^  ^+...+  2  ^  M  ;M,jj^,..,j^  positive 

integers}.  Then,  the  expectations  E^{x(d,k},  for  k  ^  0,  can  be  as  closely 

approximated  as  desired,  by  selecting  d[A  values  in  V,  in  conjuction  with  the  corresponding 

recursive  expressions.  We  do  not  include  the  latter  expressions  in  this  paper,  due  to 

lack  of  space.  The  Interested  reader  may  seek  reference  [9].  Using  the  above  methodology, 

u  •£.  u 

for  M=6  in  the  subset  V,  we  computed  the  values  of  the  bounds  H  ,  H  ,  S  ,  S  in  (B.23), 


which  are  included  in  table  B.5  below.  The  bounds  in  table  4  were  computed  from  the 
values  in  the  latter  table,  in  conjuction  with  the  expressions  in  (B.9). 


b.io 


\ 

/ 

s“ 

0.01 

1.9984 

2.0577 

0.06287 

0.06446 

0.05 

1.9244 

2.2128 

0.374114 

0.425699 

0.10 

1.96223 

2.5658 

1.00807 

1.30892 

0.15 

2.17837 

3.22266 

2.32733 

3.4354 

0.20 

2.78712 

4.6554 

6.16544 

10.2939 

0.25 

5.1692 

9.7539 

31.1957 

58.7521 

0.27 

19.0408 

122.667 

242.601 

0.28 

18.4793 

37.9514 

492.076 

997.815 

Table  B«5 


Bounds  on  S  and  H  for  the  LFS  Algorithm 
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